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Abstract 

It is well-known that the main difficulties associated with the study of 
initial-boundary value problems for linear PDEs is given by the presence 
of unknown boundary values in any method of solution. To deal effi- 
ciently with this difficulty, we have recently proposed two alternative (but 
interrelated) methods in Fourier space: the Analitycity approach and the 
Elimination by Restriction approach. In this work we present the Analyt- 
icity approach and we illustrate its power in studying the well-posedness 
of initial - boundary value problems for second and third order evolution- 
ary PDEs, and in constructing their solution. We also show the connec- 
tion between the Analyticity approach and the Elimination by Restriction 
approach in the study of the Dirichelet and Neumann problems for the 
Schrodinger equation in the n- dimensional quadrant. 



1 Introduction 

It is well-known that the main difficulties associated with Initial-Boundary Value 
(IBV) problems for linear PDEs of the type 



£(V,dii(x,t) = /(x,i), «(x,0)=« (x), xeVcK n , t>0, (1) 
at 

where V = ••, gf - ); £ is a constant coefficients partial differential operator, 
u{x.,t) is the unknown field, /(x, t) is a given forcing and ito(x) is the given 
initial condition, with Dirichelet, or Neumann, or Robin, or mixed, or periodic 
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boundary data on dV, is given by the presence of unknown Boundary Values 
(BVs) in any method of solution. To deal efficiently with this difficulty, we have 
recently proposed two alternative (but interrelated) methods in Fourier space: 
the Analyticity approach and the Elimination by Restriction approach. 

The first step, common to both methods, consists in rewriting the PDE (1), 
defined in a space-time domain T>, in the corresponding Fourier space, using the 
Green's formula. The PDE in Fourier space takes the form of a linear relation 
among the Fourier Transforms (FTs) of the solution, of the initial condition 
and of a set of BVs, only a subset of which is given a priori. This relation is 
always supplemented by strong analyticity requirements on all the FTs involved, 
consequence of the geometric properties of the space-time domain V. 

The second step is where the two methods separate; once the problem is 
formulated in Fourier space, we propose the following two alternative strategies. 

i) The Analyticity approach, which consists in using systematically the analyt- 
icity properties of all the FTs involved in the above relation, to derive a system 
of linear equations which allows one to express the unknown BVs in terms of 
the known ones, and therefore to solve the problem. 

ii) The Elimination by Restriction (EbR) approach, which consists, instead, in 
applying to the above linear relation in Fourier space a suitable annihilation 
operator, which eliminates all the unknown BVs, generating a new transform, 
well-suited to the specific IBV problem under scrutiny. The inversion of this 
new transform (if it exists) leads to the solution. 

The Analyticity approach is inspired by Fokas' recent discovery of the global 
relation, obtained first within the x — t transform approach [1] and more recently 
using differential forms [2]. The use of the global relation to study the well- 
posedness and solve IBV problems is illustrated, for instance, in [3], [4], [5]. 
In [5], in particular, general results on the well-posedness of IBV problems for 
dispersive 1 + 1 dimensional equations of arbitrary order are announced. 

Our main contribution to the method consists, after formulating the IBV 
problem in Fourier space using Green's formula, in imposing systematically the 
analyticity properties of all the Fourier transforms involved in the problem, 
to derive a cascade of analyticity constraints which allow one to express the 
unknown BVs in terms of the known ones, and therefore to solve the problem. 
In particular, Fokas' global relation appears, in the methodology we propose, as 
a "zero residue condition" for the FT of the solution. 

The Analyticity approach in the form we propose is very elementary and, 
above all, has the great conceptual advantage to originate from a single guiding 
principle: the analyticity of all the Fourier transforms involved in the problem. 
It is the type of approach that can be easily taught in elementary University 
courses, combining nicely standard PDE theory tools, like the Green's formula 
and the Fourier transform, with elementary notions in Complex Functions the- 
ory. 

The essential aspects of the Analyticity approach were first presented by 
the authors at the Workshop "Boundary value problems" in Cambridge, De- 
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cember 2001, inside the programme: "Integrable Systems". The method in its 
final form is presented for the first time in this paper, illustrated on the study 
of IBV problems of various type (Dirichelet, Neumann, mixed, periodic) for 
some second and third order classical PDEs of the Mathematical Physics: the 
Schrodingcr, the heat and the linear Korteweg-de Vries equations. Also its con- 
nections with the EbR approach are illustrated in this work, on the particular 
example of the Schrodingcr equation in the n-dimcnsional quadrant. A good 
account of the EbR approach is given instead in [6]. A different approach, valid 
for semicompact domains, has been recently presented in [7]. A general review 
of the basic spectral methods of solution of IBV problems for linear and soliton 
PDEs is presented in [8]. 

§2 is devoted to the presentation of the Analyticity approach, while §3 is 
dedicated to its application to the solution of some IBV problems for second 
and third order evolutionary PDEs in 1 + 1 and in n+ 1 dimensions. §4 is finally 
devoted to the study of the connections between the Analyticity approach and 
the EbR approach. 

2 The Analyticity Approach 

2.1 The Fourier transform and its analiticity properties 

The natural FT associated with the space - time domain V = V ® (0, oo) (in 
short: FT-p) is defined by 

F(k,q)= [ dxdte- i ( k - x+9t >F(x,t) (2) 

for any smooth function F(x, t), (x, t) G V, assuming that F(x, t) — > 0, t — > oo 
fast enough; where k = (fei, ..,k n ) € 1Z n , g £ K and k • x = Y^j=i hj x j- Its 
inverse: 

F(x,i) Xl >(x,i) = J^ +i J^e^^F^q) (3) 

reconstructs F(pc,t) in T> and zero outside, where Xvfat) is the characteristic 
function of the domain T>: xx>(x,i) = 1, (x, t) G T>, xx>(x, t) = 0, (x, t) V 
(therefore: xx>(x, t) — xv( x )H(t), where H{t) is the usual Heaviside (step) 
function) . 

If the space domain is the whole space: V — lZ n , the FTx> (2) is defined in 
A = lZ n ® T q , where 2 q is the closure of the lower half q-plane T q , analytic in 
q G X q , Vk G TZ n and exhibits a proper asymptotic behaviour for large q in the 
analyticity region. If the space domain V is compact, the FT-p acquires strong 
analyticity properties in all the Fourier variables: it is defined in A = C n <S> T q , 
analytic in q G l g , Vk G C n , entire in every complex kj, j = 1, .., n \/q € T q 
and exhibits a proper asymptotic behaviour, for large (k, q), in the analyticity 
regions. If the space domain is semi - compact, then the analyticity in the 
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Fourier variables kj, j = 1, .., n is limited to open regions of the complex plane, 
depending on the geometric properties of the domain V. 
We are therefore led to the following definition: 

Definition of admissibility. Given a space-time domain V, a function of 
(k, q) is an admissible Fourier transform for the domain V (an admissible 
FTt>) iff it possesses the analyticity properties and the asymptotic behaviour 
corresponding to that domain. 

2.2 The IBV problem in Fourier space 

We find it convenient to rewrite the IBV problem (1) in Fourier space. This goal 
is conveniently achieved using the well - known Green's formula (identity): 

bCa — aCb = div J(x, t), (4) 
and its integral consequence, the celebrated Green's integral identity: 

/ (bCa — aCb)dxdt = / J(x, t) ■ vda, (5) 

JV JdV 

obtained by integrating (4) over the domain D and by using the divergence 
theorem. In equation (4), C is the formal adjoint of C: C = C(—\J, — Jj), J(x, t) 
is an (n + l)-dimcnsional vector field, div is the (n + l)-dimensional divergence 
operator and a(x, t) and &(x, t) are arbitrary functions. In equation (5), da is 
the hypersurface element of the boundary and v is its outward unit normal. We 
remark that, given C, its formal adjoint C and two arbitrary functions a and 
b, an (n + l)-dimensional vector field J(x, t) satisfying the Green's formula (4) 
always exists and can be algorithmically found to be a linear expression of a, 6 
and their partial derivatives of order up to N — 1, if C is of order N. 

The arbitrariness of a and b allows one to extract from (4) and (5) several 
important informations on the IBV problem; with the particular choice 

a = u(x,t), b = e- i Q« x+qt )/C(ik,iq), (6) 

where C(ik,iq) is the eigenvalue of the operator £, corresponding to the eigen- 
function e 4 ( kx +9*) j the vector field J takes the following form: 
J = e ~''( k ' x +'? t ) J'(x, t; k, q)/C(ik, iq) and the Green's integral identity (5) gives 
the FTjj of the solution in terms of the FTp's (or, maybe, of generalized FT's) 
of the forcing and of all the IBVs: 

U ^'^" £( jk ,zg) -^ik,^)' 

(7) 



If the PDE has the following evolutionary form: 

dt' ~ dt 



Av,|) = |r-/C(v) (8) 
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then 

-rv \ /(k,g)+a (k)+g(k, g ) Af(k,g) 

u(k ' ?) = Z(fl^ =: Z^' (k ' 9)e ' 4 (9) 

and the linear relation (9) makes clear how the different contributions coming 
from the equation (the denominator C), from the forcing /, from the initial 
condition ii and from the set of boundary values B separate in Fourier space. 

Its inverse transform (3) gives the corresponding Fourier representation 
of the solution: 

U(x, t) = »(*,*)»<*.*) = / ^e-*"«>^, (*,<) e K-» 

(10) 

Two sources of problems arise at a first glance of equation (7): 

i) the RHS of the equation depends on known and unknown BVs; 

ii) apparently the RHS of the equation is not an admissible FTx> ■ 

It is very satisfactory that the analyticity constraints which make the RHS 
of (7) an admissible FTx> provide also a number of relations among the IBVs 
which are sufficient to express the unknown BVs in terms of known boundary 
data. 



2.3 The analiticity constraints and their resolution 

In general, C(ik, iq), the denominator of equation (7), is an entire and, most fre- 
quently, polynomial function of all its complex variables. Let S be the manifold 
in which this entire function is zero: 

S = {(k, q) e C n+1 : £(ik, iq) = 0}. (11) 

Then the RHS of equation (7) provides an admissible FTx> of the solution of the 
IBV problem under investigation if the numerator N (k, q) of u in (7) satisfies in 
AC\S, hereafter called the singularity manifold (SM) of the IBV problem, 
the following Zero Residue Condition (ZRC): 

A r (k,g) = 0, (k,q)eADS. (12) 

If the singularity manifold AdS contains the real axis (which is usually a part 
of the boundary of A) and if this singularity is not already taken care of by the 
ZRC (12), then we must also proceed to the Denominator Regularization 
(DR): 

C(ik,iq) -> Cregiik, iq), (13) 

which consists in moving a bit the singularity off the real axis, outside the 
domain A. 

The ZRC plus the DR constitute the main set of Analyticity Constraints 

(ACs) that must be imposed to the RHS of (7) in order to obtain an admissible 
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FTt> of the solution of the IBV problem under investigation. The ZRC (12) 
provides a (linear) relation among the FTs of the forcing, of the initial condition 
and of all the BVs. The analyticity properties of all these FTs generate, through 
the admissibility argument, a cascade of further analyticity constraints, untill 
all these conditions are finally met. This procedure defines, in principle, a set 
of relations (a system of equations) among the IBVs. Therefore: 

a) The unique solvability of such a system, together with the admis- 
sibility of the obtained solution, are equivalent to the study of the 
unique solvability of all the IBV problems associated with (1). 

b) By solving this system for a set of BVs in terms of the remaining 
ones, one expresses all quantities in terms of known data and, from 
equation (10), one obtains the Fourier representation of the solution. 

In most of the examples considered in this paper, this system of equations 
is algebraic, with entire coefficients. Therefore, if M is the squared matrix of 
the coefficients of the unknown BVs, the admissibility argument imposes that 
the countable number of zeroes of det M: 

det M(q m ) = (14) 

lie outside the analyticity domain of an admissible FTv- 

q m il q , me TV. (15) 

It turns out that the set (14) coincides with the spectrum arising in the eigen- 
function expansion approach [9] and coincides also with the restricted domain 
in which the EbR method works. These deep connections justify for (14) the 
name of spectrum of the IBV problem. 

The admissibility argument imposes also that the constructed solution of the 
system exhibit the proper asymptotic behaviour in the analyticity domain. It is 
actually convenient to impose first this asymptotic admissibility, the easiest to 
be checked, which enables one to disregard without effort all the IBV problems 
ill-posed because inconpatible with asymptotics. 

2.4 General remarks 

Remark 1. Analyticity vs Causality. It is well - known that there are 
definite connections between the analyticity properties of the FT of the solution 
of evolution equations and the causality principle. In our general setting it is 
straightforward to show that: 

The analyticity properties of the FT?, of the solution of the IBV prob- 
lem (1) imply the causality principle. 

Indeed, using the convolution theorem, the inverse FT (10) of the RHS of equa- 
tion (7) (in which all the analyticity constraints have been preliminary imposed) 
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is equivalent to the following Green's representation of the solution: 



u(x,t) = J dt 1 j ' dx.'GRF(x-x';t-t')Af{x' ' ,t'), (x,i)eD, (16) 

V 

where 7V(x, t)xv(x, t) is the inverse FT (3) of 7v(k, q) and Grf is the celebrated 
retarded - fundamental Green's function of the operator C: 

G RF {*,t = / - ; 17 

(2tt)" £ re9 (zk,2g) 

which satisfies the important property: Gnp{x.,t) = 0, i < 0, due to the reg- 
ularization of £(ik,iq). Equation (16) is the usual way in which the causality 
principle becomes transparent. 

Remark 2. Regularization and Fourier representation. As we have al- 
ready written, if the zeroes of the denominator on the real axis are all cured by 
the ZRC, no regularization is needed. On the other hand, some regularization 
must be introduced also in this case, in the calculation the Fourier represen- 
tation (10), before splitting N in the sum of the different contributions (each 
one singular on the real axis) in (9) coming from the forcing, from the initial 
condition and from the BVs. The most convenient regularization is obviously 
that in (13) and leads to the following Fourier representation: 

c rfkd ge '' kx +^> go(k) i c dkdqe< kx+ ^ B(k,q) / vric ™tl ^ ' 

JK"+1 (2tt)"+i C rcg (ik,iq) + JK"+1 (2tt)" + i C rcg {ik,iq) ' l^, t /-C 



3 Illustrative Examples 

In this section we apply the Analyticity approach to the following classical 
equations of the Mathematical Physics, the Schrodinger, the heat and the linear 
Korteweg-de Vries (KdV) equations: 

du O^u 

-di~ a d^ = ^ a = hl x ^ v > t>0 > ( 19 ) 
du d 3 u 

dt~ V dx^ = L V = ±1, X £ V > 1 > °' (20) 
prototype examples respectively of second and third order evolutionary PDEs 
and basic universal models for the description of dispersive and diffusive phe- 
nomena, where the space domain V is either the segment (0, L) or the semiline 
(0, oo). Hereafter the BVs will be indicated by 

v { o j \t):=^(x,t)\ x=0 , vP(t):=^(x,t)\ x=L , jeM (21) 
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and their Fourier transforms by Vq\q), v^\q)'. 

oo 

« W) (g) := / dte-^v^it), v^(q) := / die"^^*). (22) 



We also apply the method to the study of IBV problems for the multimen- 
sional analogue of equation (19), for a = i: 

du d 2 

— -iAu = /, xeV,t>0, A:=w=Ep ( 2 3) 

j=i 3 

in the n-dimcnsional Quadrant 

y = {x: >0, j = l,..,n}. (24) 
The corresponding BVs will be indicated by: 

v$(xj,t) = u(x,t)| x , =0) ^(xj-.t) - ^L(x,t)| x . =0 (25) 
and their FTs by: 

4™ ) (k j , g ) = /^/rfx je -^^+9*)^™ ) (x j ,i), m = 0,l. (26) 



In equations (25)-(26) Xj = (xi, .., Xj, .., x„) G 7£ n_1 , kj = (fci, .., fcj, .., fc n ) G 

7£ n_1 , J y c?x^ = J Ll dx' x • -(/ ij dajj) • • J Q L " e^4 and the superscript " indicates 
that the quantity underneath is removed. 

The application of the Analyticity approach to higher order problems and 
to other relevant examples will be presented in [8] . 

3.1 The second order PDEs (19) 

In this case, equations (4) and (6) imply: 



(27) 



C(ik, iq) = i(q — iak 2 ). 
In addition, if V is the segment (0, L), equation (9) yields 

N(k, q) = f(k, q) + fio(fc) - a {[v^ (q) + ikv^ (q)} - e~ ikL [v™ {q ) + lk vf (g)]) 



(28) 
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and the Fourier representation (18) of the solution takes the following form: 

«(*• *) - /*» (S 6 ^^^ ^fi^lo + J* ^-«* , %(t)+ 

Jac<°> f-e^-^ 2 '^ 1 ^^ 2 ) +^4 0) (*« fc2 )]+ (29) 
/ac<"> |^e lfc ( a: - L )- Qfe2t [4 1) (*afc 2 ) +^4 0) (*a fc2 )], & G (0,i), * > 0, 

where JC^ and /C^ are respectively the first and third quadrant of the complex 
k - plane, K, 1 ^ = p^Km , m = 0,1, where qtl is the 7r/4 rotation operator: 

p* : fc — > e^fc and 9/C is the counterclockwise oriented boundary of /C (see 
Figs la,b). 

The corresponding expressions for the semiline or for the infinite line cases, 
with rapidly decreasing conditions at oo, follow immediately from the ones 
above, setting vf ^ = = in the semiline case, or setting = = 
vf^ = = in the infinite line case. 

It is instructive to first apply the Analyticity approach to the simplest case 
in which the space domain is the whole space, with rapidly decreasing BVs at 
x = ±oo. 

3.1.1 The whole line V = (-00,00) 

Equation (28b) reduces to 

ti(k,q) = f{k,q)+u {k) (30) 

and the admissibility argument imposes that N{k,q)/C be defined in (k,q) E 
A = 1Z ®T q and be analytic in q E I q , V/c G 1Z. If a = 1, the denominator is 
singular for q — ik 2 , k G 1Z, outside the definition domain, and no rcgularization 
is needed. If, instead, a = i, the denominator is singular for q = —k 2 < 0, on 
the real negative axis, at the boundary of the analyticity domain, and the only 
analyticity constraint to be fulfilled is the Denominator Regularization (13): 

C(ik, iq) = i(q — iak 2 ) — > C reg (ik, iq) = i(q — iak 2 — iO). (31) 

The regularization (31) is sufficient to make the RHS of (28a) an admissi- 
ble FTt>, from which we recover the well-known Fourier representation of the 
solution of equations (19): 

f dgdfce'( fa +*> f(k,g) f dk , 
U{X,t) -Jv (2n)H q-iak 2 -i0 + ] n 2^ u (k) , x E K, t > 0. 

(32) 

We now proceed considering semi - compact and compact domains. 
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3.1.2 The semiline V = (0, oo) 



In this case: 




q e Iq, Vfcelfe and be analytic in k e V? £ I g . Therefore the singularity 
manifolds AP\S^ a \ corresponding to a = i, 1, are parametrizable either in terms 
of k or in terms of q in the following way: 



ADS {a) ={q = tak 2 , k e IC ( a) } = {k = k { a \q), ir < arg q < 2tt}, (34) 



If a = 1, there is no singularity on the real axis and no regularization is 
needed. If a = i, there are two singularities for k G 1Z; that corresponding to 
k < is cured by the ZRC (12), while that corresponding to k > is cured 
instead by the regularization (31). 

The ZRC (12) is conveniently parametrized in terms of q in the following 
way: 



fi(kk a) (q),q) = n^ a \q),q)+M^ a \q))-a[4 1 \q)+ i k^\q)v^(q)} = 0, 



for 7r < arg q < 2ir. It is one equation involving 4 FT' s which are therefore 
dependent. If we are interested in solving the Dirichclet and Neumann problems, 
we use this ZRC to express the unknown BVs in terms of the known ones: 



Dirichelet problem : av^^q) = f(kg a \q),q) + u (kQ a \q)) — aikQ a \q)v^\q) 
Neumann problem : iakQ a \q)v^ (q) = f(k^(q),q) + uo(kQ a \q)) —ai>^\q), 



for 7r < arg q < 2tt. It is easy to see from (37) that the unknown BVs define 
admissible FTs which, inserted in (29), give the wanted solution of the Dirichelet 
and Neumann problems. 

We remark that the ZRC (36) could also be solved for m (using now, for 
convenience, the variable A:): 



u (k) = -f(k,iak 2 ) +a[v^\iak 2 ) + ikv^ {iak 2 % k e JC ( Q a) (38) 



but, in this case, the solution would not be, in general, an admissible FT, since 
the RHS of (38) cannot be extended to the rest of the lower half k - plane. 
Even in the special case in which the forcing and the assigned BVs were on a 
compact support in t, corresponding to entire FTs, the solution uo(k) would not 



where 




(35) 



(36) 
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be admissible, because it would not possess, in general, the proper asymptotics. 
This means that the (unphysical) problem in which we assign arbitrarely u and 
its space derivative at x = cannot be treated by this method, unless the above 
BVs are suitably constrained. 

3.1.3 The segment V = (0, L) 

Now admissibility implies that Af/C be defined in A = C <£) T q , be analytic in 
q € T q , Vfc G C and be analytic in k G C, Vq G I q , with proper asymptotics 
for large |fc and/or \q\ in the analyticity regions. Therefore the singularity 
manifolds on which the ZRC (12) is defined are now the unions of two sectors: 



A n = U m=0 {g = iak 2 , k G /Cm } = 
U m=0 {fc - *&\q) = (-) m k { a \q), tt < arg q < 2n}. 

k plane k plane 



Fig.la The SM A n S (a) (a = i) Fig.lb The SM A n S (a ^ (a = 1) 

Both singularities on the real axis are cured by the ZRC and no regularization 
is needed. The regularization (31), however, is still introduced, according to the 
Remark 2 of §2.4, in computing the Fourier representation (65) of the solution. 
The ZRC (12), conveniently parametrized using q, consists of the following 
system of two linear algebraic equations: 

Af(k$ {q),q)=0, m = 0, 1 tt < arg q<2ir (40) 

containing four BVs. Therefore we expect to be allowed to assign arbitrarely two 
out of four BVs. To establish which pairs of BVs can be assigned arbitrarily, one 
should impose that the corresponding solutions of the algebraic system define 
admissible FTs; i.e., the following two conditions must be satisfied, 
i) The system must be uniquely solvable for the unknown pair of BVs in its 
definition domain. More precisely, indicating by M the 2x2 matrix of the 
coefficients of the unknown BVs, the admissibility condition imposes that the 
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countable set {qj}jej^ of zeroes of det M, the spectrum of the IBV problem, lie 
outside the analyticity domain: 

q^l q , j&Af. (41) 

ii) The solution of the system must define admissible Fourier Transforms; in 
particular, it must exhibit the proper asymptotics in the analyticity domain. 

Studying first the asymptotics of (40), one infers without any effort which 
pairs of BVs cannot be assigned arbitrarily. The asymptotics of (40) imply 
immediately that the following expressions: 

f{k<?\q), q) + u (k { o a \q)) - ^ (q) + ik™ (<f)t>< 0) (q)} , 
e ik o a) M L [f(-ki a \q),q) + u (-k ( o a) m + rffe) - 4 Q) (9)4 0) (<7)] 

are exponentially small for q ~ cxo in n < arg q < 2n. Since the asymptotic 
scries of the admissible FTs appearing in the LHS of equations (42) are inverse 
power series of <p , equations (42) impose severe constraints on the involved 
functions, implying that: 

asymptotic admissibility is compatible with assigning at x — any BV between 
(v^ , Vq 1 ^) and, at x = L, any BV between (v^\ v^). It is not compatible 
instead with assigning arbitrarily the pairs (v^^v^) or (v^\v^). 

To complete our analysis, we must check if the spectrum associated with 
the IBV problems compatible with the asymptotics lie outside the definition 
domain. The analysis is straightforward and produces the following results. 
Proposition (the spectrum). Assigning arbitrarely (v^\v^) (the Dirichelet 
problem) or (vq~\v^) (the Neumann problem), the spectrum is characterized 
by the equation sin(fcL) = k m = ^p, m G Z and is given by the 

negative eigenvalues {qm}neMi q-m = — k^ = ~(^) 2 , rn G Af, if a — i, and 
by the purely imaginary eigenvalues q m = ikf n = i(^) 2 , m G Af, if a = 1. 

Assigning instead (v^\v^) or {vq~\v^) (the mixed problems), the spectrum 
is characterized by the equation cos(fcL) = k m = |(2m +1), m G Z 
and is given by the negative eigenvalues {q m }meN, qm = — &m = — (-|) 2 (2m + 
l) 2 , m G Af, if a = i, and by the purely imaginary eigenvalues q m = ik^ = 
i(f ) 2 (2m+ l) 2 , m G Af, if a = 1. 

For a = 1 the spectrum lies outside the analyticity region and the solutions of 
the algebraic system (40) define directly admissible FTs; if a = i the solutions of 
the algebraic system (40) define admissible FTs after moving these singularities 
a bit off the real q - axis, outside the definition domain (again a regularization!). 
We conclude that all the IBV problems compatible with admissible asymptotics 
turn out to be well-posed: 

IBV problems for the Schrodinger and heat equations (19) are well- 
posed assigning at x = any BV among (vq > \vq''*) and at x = L any 
BV among (vf \v^). 
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It is interesting to remark that, if one insisted, instead, in solving an IBV 
problem in which the BVs (vQ~\q), (q)) are assigned, the correspondig al- 
gebraic system would be always uniquely solvable (no point spectrum would 
arise), but the solution would exhibit an exponential blow up at q ~ oo in the 
analyticity region, that cannot be accepted. This undesired blow up could be 
cured if the assigned BVs were related by the (additional) analyticity constraint: 

f(k { a \q),q)+u (k { a) (q)) - a[v^ (q) + ik™ (g)t>< 0) (q)} = 0, (43) 
implying the following admissible solutions of the algebraic system (40): 

v£\q) = -ik^\q)vf\q) = -^|^[/(/>) (?),</) + /("^ («),?)+ (44 ) 
tio(^ a) («))+«o(-^ a) («))-2cw)( 1) (g)]. 

The additional analyticity contraint (43) is not surprising at all, since it is 
nothing but the ZRC of the scmilinc problem. Similarly, assigning the right 
boundary conditions {vf \ v^), the unknowns and would exhibit again 
an exponential blow up which cannot be accepted; an admissible asymptotics 
would be guaranteed now by the (additional) analyticity constraint: 

e-< a) ^ L [f(-^ a \q),q)+M-^ a \Q))]+^ = 0, (45) 

implying the following solution of the algebraic system: 



vg\q) = iki a \q)vi°\q) = ^-[f(ki a \ q ),q)+f(-k^(q),q)- 
Mk^\^)+M-ki a \q)) + 2ae-< > ^ L 4\q)}. 

3.1.4 The periodic problem 



(46) 



~(o) -(o) 
— V, 



If we assume L-periodicity of u and u x , then Vq 1 ^ = = : Wq° - ■ . 
t)(°) and the algebraic system (40) consists now of two equations for two BVs, 
which have to be treated therefore as unknowns. The solutions of this system 
read: 



~(1) /„-> _ 1 ( /(fc^ ) (g),g)+flo(fcf,° ) ('?)) , /(-4° ) (9)^)+"o(- fc o° ) (g)) 

(g),g)+fio(-fr (a) ' 



f,(0)(V| = 1 f /(C(g).9)+""(C J ('?)) _ f(.-K c "(q),q)+u a (-k^>(q)) 



(47) 



They satisfy asymptotic admissibility and the spectrum, characterized by the 
equation 1 — e ±lkL — 0, (=> k n — ^n, n 6 Z), is given by q n = —k\ = 
-(^) 2 n 2 , n e Af, for a = i, and by q n = ik 2 n = i(^-) 2 n 2 , n 6 TV, for a = 1; 
therefore the usual regularization is needed again in the Schrodinger case. We 
conclude that 
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the periodic problem for equations (19), in which one imposes the 
L- periodicity of u and u x , is well-posed and no BV can be assigned 
arbitrarely. 

Remark We remark that the Fourier transforms of the unknown boundary 
functions exhibit generically a branch point at q = 0, due to the well- known 
slow decay as t — > oo of the solutions of the dispersive evolution equation under 
investigation. 

The above procedure generalizes with no difficulties to higher order prob- 
lems. In the following we concentrate on a third order problem only. 



3.2 The linear KdV equation 

In this section we investigate IBV problems for 3 rd order operators, illustrating 
the method on the simplest possible example (20). 

Since the group velocity v g = 3i]k 2 of the associated wave packet is positive 
(negative ) for rj positive (negative), we have the following expectations. In the 
scmiline case, one should be able to assign at x — more BVs for positive r\ 
than for negative r\. In the segment case, for r\ positive one can assign arbi- 
trarely more BVs at x — than at x = L (and viceversa for r\ negative). The 
precise indication of "how many" and "which" BVs can be assigned in order to 
have a well-posed IBV problem follows again in a straightforward way from the 
Analyticity approach. 

Equations (4) and (6) imply: 



(48) 



C = -C, J=(a6,-#^-f f + f£o]), 
£(ik, iq) = i(q + rjk 3 ). 

In addition, if V is the segment (0, L), equation (9) yields 

Af(k,q) - f(k,q) + U (k) - V {[vP(q)+ik4 1 \q) - Pv^(q)}- (49) 
e- ikL [v?\q)+ikv£\q)-k 2 vW(q)}) 

and the Fourier representation (18) of the solution takes the following form: 

«(*.*) = /*■ ^ i(kX+qt) ^Fo + k ge^-^fioflfc)- 
v (f^ n) ge*^-"* 3 *)^ 2 ^-^*; 3 ) + ikv { o\-r]k 3 ) - k 2 v { ° } (- v k 3 )}+ 
j ^e i ^ x -° L ^ k3t ^v^\-r ] k 3 )+ikv { L\-r 1 k 3 ) - Pv^ (-r)k 3 )}), x e (0,L), t > 0, 

where 7o " = dJC { -\ ^ = dK^ U 8K { ~\ 7 + = 9/c[ +) U dK, { +\ 7 + = dK, ( +\ 

Ktf = {k : J(2m + 1) < arg k < ^(2m + 2)}, /C« = m = 0, 1, 2, 

(51) 
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and /3,r is the involution p n : k — > —k. 
3.2.1 The segment V = (0, L) 

Now A/"/£ must be defined in A = C ®I q , analytic in 5 e 1,, Vfc G C and 
analytic in k e C, Vg £l 9 , with proper asymptotics for large |fc| and/or g| in the 
analyticity regions. Therefore the singularity manifolds AnS^ v \ corresponding 
to r\ = ±1, are given by (see Figs 2a, b): 

.AHSW = Ul l=0 {q = - V k 3 , k £ = U 2 m=Q {k = k%\q), tt < arg g < 2n} 

(52) 

where km\q) = —r)p m q^ and p TO , m = 0, 1, 2 are the 3 roots of unity: 

p m = e^ m , m = 0,l,2. (53) 



k plane k plane 




Fig.2a The SM AnS™ (rj = -1) Fig.2b The SM A D <S^ (77 = 1) 

The ZRC (12) consists of the following three equations: 

Af(fcW(g),g) = 0, m = 0,l,2, tt < ar ff g < 2tt. (54) 

For q E 1Z there is one singularity on the real k - axis, which is cured by one 
of the three equations (54) and no denominator regularization is then needed. 
The regularization (31), however, is still introduced, according to the Remark 
2 of §2.2.2, in writing the Fourier representation (50) of the solution. 

The 3 algebraic equations (54) contain 6 BVs; therefore we expect to be 
allowed to assign independently only 3 BVs. As before, a quick asymptotic 
extimate selects the sets of 3 BVs which can be assigned independently, com- 
patibly with asymptotic admissibility. The asymptotics of equations (54) imply 
that the following expressions, respectively, for r\ = — 1: 

e^ L [f(qi,q)+Mqi)} - [«£ 2 > (g) + iqi (g) - qiv f (g)], 

Hpii^q) + u {piq^) + [4 2) («) + l / , i9^o 1) (9) ~P2g^o 0) (?)], ( 55 ) 
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and for rj = 1: 

f(-Q*>Q) +uo(-q*) - [i> ( o\q) - iqh^ } (q) -qiv^ 0) {q)], 

e-^^[f(- Piq i, q ) + u (-piq*)} + [v { L ] {q) - ipiq^Hq) - m^\ q % (56) 

e- i "^ iL [f(- P 2qKq) + «o(-p 2 ^)] + [v?\q) - im^£' '(«) - Pi^f (?)], 

are exponentially small for q ~ oo in 7r < arg g < 27r. Therefore, reasoning as 
before, we see that: 

i) for 77 = — 1, a necessary and sufficient condition to obtain FTs with admissible 

asymptotics is to assign at x = any BV among v^\v^\v^ (consequence of 

equations (55b, c)) and, at x = L, any two BVs among vf \v^\v^ (conse- 
quence of equation (55a)); 

ii) for 1] = 1, a necessary and sufficient condition to obtain FTs with admissible 
asymptotics is to assign at x — any two BVs among , , (consequence 
of equation (56a)) and, at x = L, any BV among vf\v^\ (consequence of 
equations (56a, b)). 

Again, to complete our investigation, we must check if the spectrum associ- 
ated with the above IBV problems selected by the asymptotic admissibility lie 
entirely outside the analyticity domain of an admissible FT. It is easy to prove 
that it is indeed the case. 

Proposition (the spectrum of the IBV problem) Consider any IBV prob- 
lem on the segment for equation (20) compatible with the asymptotic admissibil- 
ity established above; i.e., in which, for n = —I, one assigns arbitrarely at x = 

DI , (0) (1) (2) , , DTA , T (0) (1) (2) 

any B V among Vq ,Vq ,Vq and any two BVs at x = L among v L , v L , v L , 
and in which, for n = 1, one assigns arbitrarely at x = any two BVs among 
i>Q°\i>Q 1 ' ) ,t>Q 2 ' ) and any BV at x = L among vf\ v£\ . For rj = —1, let Vq 1 ^ 
be the given BV at x = and be the unknown BV at x = L while, for 

rj = 1, let Vq be the unknown BV at x = and v ™ be the given BV at x = L . 
Then the corresponding spectrum is characterized by the following equation: 

A("( m ""))(fc) = 0, (57) 

where: 

A^(fc) := e- lkL + p{e- plikL + p{e- p2ikL . (58) 

The proof is tedious but straightforward and makes essential use of the well- 
known algebra of the roots of unity, which implies also that all the above 
IBV problems lead only to three (similar) purely imaginary discrete spectra 
{kn^}neAf, characterized by the three equations A^\k) =0, j = 0,1,2. More 
precisely: 1) the spectrum characterized by equation A^°^(fc) = is given by: 

k^ = -i{^/L),n&N+: {d 0) W + : 4^ = - cos(f m 
=> d 0) ~^(2n-l), n>l. 
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2) The spectrum characterized by equation A^(() = is: 

=► d X) =0, j£(n-§), n>2. 

3) The spectrum characterized by equation A^(() = is: 

fci 2) = -i^/L), neN: {d 2) W : = cos(^ 2) - f ), 

C f =0, Ci 2) ~^(n-J), n>l. 

We conclude that all the three discrete spectra 



(60) 



(61) 



{AeM, qP=W 3 =i(^-) , .7=0,1,2, (62) 



associated with the above IBV problems lie on the positive imaginary axis of the 
complex q plane, outside the analyticity domain of an admissible FT. Therefore: 

IBV problems for equation (20) on the segment (0, L) are well-posed 
iff: 

i) for 7] = —1, one assigns at x = any BV among v^, and at 

x = L any two BVs among vf\v^\vf^; 

ii) for rj = 1, one assigns at x — any two BVs among Vq~\ and 
any BV at x = L among 



3.2.2 The periodic problem 

If we assume L-periodicity of u, u x and u xx , then = j = 0,1,2, the 
algebraic system (54) consists now of three equations for three BVs, which have 
to be treated then as unknowns. The solution of this system satisfy asymptotic 
admissibility and the spectrum, characterized by the equations 1 — e~ Vjk — 
0, j = 0,1,2 (=> k n = ^-pj 1 ™., n e Z, j = 0,1,2), is given by the real 
numbers q n = —rjk^ = — ry(^-) 3 n 3 , n € Z and must be regularized in the usual 
way. We conclude that: 

the periodic problem for the linear KdV equation (20), in which one 
imposes L-periodicity to u, u x and u xx , is well-posed and no BV can 
be assigned. 



3.2.3 The semiline V = (0, oo) 

Taking the limit L — > oo of the results of §3.2.1 we immediately obtain the results 
on the semiline. In this case, the singularity manifolds are the restrictions of the 
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above ones to the lower half k plane. No spectrum arises and the asymptotic 
admissibility implies that: 

IBV problems for equation (20) on the semiline (0, 00) are well-posed 
iff, for T] = —1, one assigns at x = any BV among (uq°\ v o^\ w o^) an d, 
for r] = 1, one assigns at x = any two BVs among (vq*\vq~\v^). 

We remark that, in the cases treated so far, the spectra of all the IBV prob- 
lems compatible with asymptotic admissibility lie always outside the analyticity 
domain. We do not have, however, a general argument excluding the situation 
in which part of the spectrum lie inside. Therefore the complete characteriza- 
tion of the spectrum, the only part of the method in which some technicality is 
involved, seems to be unavoidable and makes it difficult to prove general results 
for operators of arbitrary order. 

The Analyticity approach applies nicely also to an arbitrary number of di- 
mensions and next section is devoted to an illustration of it. The application 
of the method to higher order problems and to other relevant examples will be 
presented in [8]. 



3.3 Multidimensional Schrodinger equation 

In this section we study the Dirichclet and Neumann problems for the Schrodinger 
equation (23) in the n-dimensional quadrant (24). Then: 



C = — ^ — iA, J = (ab, i(a\/ b — bxj a), 
£{ik,iq) = i(q + k 2 ), 



(63) 



where k 2 = k • k. Equations (7) and (63) give the following expression of the 
Fourier transform of the solution in terms of the Fourier transforms of the forcing 
and of all the IBVs: 



u(k, q) 



Af(k, g ) 
l(q+k 2 ) 



Af(k, q) := /(k, q) + fi (k) - i £ [v$(kj,q) + ik^f (kj , q)]. 

J'=l 

The Fourier representation (18) of the solution reads: 



«(x,t) = J nn+1 ^je^^J^ + J^ ^ e ^-^ o( k) + 
where dkj = dk\..dkj..dk n . 
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In view of the distinguished parity properties of the Fourier transforms in 
(64), we shall make an extensive use of the parity operators: 

n n 

A ± =n(i±a,), *g =n(i±m), (66) 



1=1 



where Uj is the involution Uj : kj — > — kj . 

In this multidimensional case, the FT of the solution is defined in A = 
iki ® ■ ■ ®^k n ® i g , analytic for q G T q , Vk G X kl ® • • ®2fc n and in kj G Ife 3 . , 
Vkj G Ifcj €5 • • ®Tkj <S> ■ • ®Ik n and Vq el,. We found it convenient to study the 
ZRC in the n different regions QJ C A n S, j — 1, .., n defined by: 

Qj := {(k,g) G C" +1 : k, G ft"" 1 , tt <arg q < 2tt, fc, = Xj-(kj,«)}, 
Xj(kj,g) —i^ + kj - kj-)' ei, j = l,..,n. 

Therefore the starting point of the analysis is the set of n equations 

Af(k,q)\ kj=X] = 0, kj G ir~\ tt < arg q < 2tt, j = 1, ..,n. (68) 



(67) 



Dirichelet problem. The parity properties in k of the BV terms imply that the 
application of the parity operator defined in (66) to the j th equation (68) 
eliminates all the v^s except vf£): 

AV>v$(kj,q) = -A^(W(k,q)\ k ^), 3 = h-,n, 

W(k, q) := /(k, q) + iu (k) + i £ Mo?( k ^<?) (69) 

and the analyticity properties of the v^s allow one to express them in terms of 
known quantities: 

«#(]£,•,</) =V^A^v£ > (\L j ,q) = -V^A { l\w(k,q)\ k]=K] ), j = l,..,n (70) 



applying the lower half plane analyticity projectors in all the k- variables (except 

J_ f dk' m 
2m J k' m - [k m - iO) ' 



^=f[V m , V m = -±[ v _^_ iny (71) 



Equations (70) summarize all the analyticity informations contained in the ZRC, 
allow one to express the unknown BVs in terms of given data and, via (65), to 
solve the Dirichelet problem. 

Neumann problem. Similar considerations can be made in solving the Neumann 
BV problem. In this case: 

^•(k^^A^fig^k,-^) = -A«\v(k,q)\ kj=Xj ), 3 = 

V(k, q) := /(k, q) + i«o(k) + £ ^(k,, q) (72) 
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and 

v${kj,q) = iVW (^-^Af\v(k, q )\ k]=Xj )) , j = !,..,„. (73) 

In this multidimensional context, for the presence of the analyticity projec- 
tors, the unknown BVs in Fourier space turn out to be nonlocal expressions 
of the given data. It is however possible to show that, due to the analyticity 
properties of the involved FTs, it is not really necessary to apply the above 
analyticity projectors to construct the unknown BVs and the solution it(x, t) in 
configuration space. The strategy to avoid unpleasent nonlocalities is outlined 
in the next section and leads to a Fourier representation of the solution already 
obtained in [6] using the EbR approach. Therefore this strategy is also the way 
to establish the connection between the Analyticity and the EbR approaches. 



4 Connections between the Analyticity and the 
EbR approaches 

Dirichelet problem 

We first remark that the unknown BVs can be constructed directly in terms of 
known data from the RHS of (69b): 

vS>(^,t) = -j nn ^e^ + rt A U) {W{] ,, q) \ kj=Xj) , t>0, Xk >0, k*j. 

(74) 

Indeed, from the analyticity properties of v^j we know that its inverse FTs (3) 
is zero outside the domain of definition in configuration space (i.e., for x k < 
0, k ^ j); this implies the formula 

/ dk 3 dqe^^ + ^[vff(k 3 ,q)-A^vff(k 3 ,q)] = 0, t > 0, Xj > 0, j = 1, .,n 

(75) 

and, through (69a), equation (74). 

Also the solution u(x, t) can be reconstructed without going through the 
nonlocalities associated with the analyticity projectors. Indeed it is possible to 
show that the following relation holds true: 

n 

^««(k jl9 ) = (A_-l)W(k, g ), (76) 

where the equivalence A(k, q) = B(k, q) means that the FTs A(k, q) and B(k, q) 
are equal under the following Fourier integral projector: 

J dkdq ^* 2 +qt M [i(k, q) - B(k, q)} = 0, (x, t) E T>. (77) 
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The equivalence (76) and equation (69) imply 

Af{k,q) = A-W{k,q) (78) 
and the following spectral representation of the solution in terms of known data: 

j = l K n-l 1 

already obtained in [6] using the EbR approach. 

The proof of (76) is based on the important fact that all the admissible Fourier 
transforms A(k, q) under consideration satisfy the equivalence 

o-ji(k, ? ) = i(k,g)| fe . =x ., j = l,..,n (80) 

and goes as follows. For n = 2, the 2 ZRCs (68) and their consequences (69) 
yield the 4 equivalence relations: 

Voi(k 2 ,q)+a 1 v^(k 1 ,q)= - oiW(k,g), o^oi^, g) + v^^i, ?)= - ^W{k,q), 
(1 - <r 2 )^ 1 1 ) (fc 2 , q)= - cti(1 - <r 2 )^(k, g), (1 - ai)*^ (Ai, g)= - a 2 (l - CTi)W(k, g) 

and their sum is exactly equation (76). To generalize this result to the case of 
an arbitrary n, consider the n ZRCs (68) and all their consequences, obtained 
applying systematically parity operators characterized by different indeces: 

Vo~j > Q s -j,q)+VjT,v§ ) Qs.i,q) = -&jW(k,q), j = l,..,n, 

(1 - CTjJwg^kj-.g) + - &i) v$(ki,q) = 

W (82) 
-<7j(l - &i)W(k,q), i^j, j = l,..,n, 



AV%>{kj,q) = -vj^Wfrq), J = l,..,n. 

The sum of all these equations with weights 1/ (J 1 ^ 1 ) (m is the number of parity 
operators appearing in the equation) yields the result (76). 
Neumann problem 

Similar considerations can be made in the case of the Neumann IBV problem. 
Now the unknown BVs are recovered via: 

fffte'*) = ~j nn ^e i{k ^ +qt) (^V(k, q )\ k]=X] ), t>0, Xk >0, k^j 

(83) 

and the spectral representation of the solution reads: 
«(x,t) = -J nn+1 j^ne^^^im + I^ ^e^-^A + u Q (k)- 

t I ^(XcW^e^-^Ag^yO^,-* 2 ), (x,t)eD; (<S4 ' 

J=lK»- 1 1 
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a formula already derived in [6] using the EbR approach. 
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